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Abstract 

We examine the possibility for precise determination of |Vts| (and |Vtd/^ts|) from the 
exclusive decays, B -> Kt+£~, B -» K*l + l~ (and B -> vr^", 5 -► We 
show that the ratio |Vt s | can be extracted experimentally with a small theoretical 
uncertainty from hadronic form-factors, if we appropriately constrain kinematical 
regions of q 2 . We also give detailed analytical and numerical results on the differential 
decay width dT(B — ► K*l + l~)/dq 2 , and the ratios of integrated branching fractions, 
B{B — » p£ + £~)/B(B — » We estimate that one can determine the ratio |Vt s | 

from those decays within theoretical accuracy of ~ 10%. 
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1 Introduction 



The determination of the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix is 
one of the most important issues of quark flavor physics. The precise determination of Vtd 
and V u b elements is of special significance, since they are closely related to the origin of CP 
violation in the Standard Model (SM). Furthermore, the accurate knowledge of these matrix 
elements can be useful in relating them to the fermion masses and also in searches for hints of 
new physics beyond the SM. For these reasons many strategies for the accurate determination 
of Vtd and V u b are under intensive investigation The main reason we are interested in 
B physics is that this area is very likely to yield information about new physics beyond 
the SM. We expect that new physics will influence experimentally measurable quantities in 
different ways. For example, most of us expect that AB = 2 transitions are more sensitive 
to new physics than decay rates. New physics may couple differently to K mesons compared 
to B mesons. Therefore, it is essential to determine the CKM matrix elements in as many 
different methods as possible. 

In the existing literature, we can find several proposals of different methods for precise 
determination of V td and/or \Vu/V ts \ % @, % [|: 

• | Vtd | can De extracted indirectly through Bd — Bd mixing. However, in B& — Bd mixing 
the large uncertainty of the hadronic matrix elements prevents us from extracting CKM 
elements with good accuracy. 

• A better extraction of |Vtd/Vf S | can be made if B s — B s is measured as well, because 
the ratio f B B B J p B B Bs can be predicted much more reliably. 

• The determination of |Vtd/Vt s | from the ratios of rates of several hadronic two-body 
B decays, such as T(B° -> K*°K°)/T(B° -> <f>K°), T(B° -> K*°K* )/T(B -> </>^*°), 
T(B + -> K*°K+) /T(B+ -> (f)K + ), and T{B + -> K*°K* + )/Y{B+ -> <f)K*+) has also 
been proposed in 0. 

• Vtd can be determined from K — > rruu, B — > ttujj and B — > pvv decays with small 
theoretical uncertainty |], 

• In H| a new method was proposed for the determination of |Vtd/Vt s | from the ratio of 
the inclusive decay distributions 

^-(B^x d eu-)AB^xj+n , 

as as 
where s is the dilepton invariant mass. 



2 



In this work, in order to determine \Vt s \ (and the ratio | V*rf/ V* s | ) we carefully examine 
another well known method from an analysis of exclusive decay distributions 

vr£+r), ^-(B^Ki+r), ^-(B^p£ + n and ^-{B -> K*£+r). 
as as as as 

It is well known that the experimental investigation and detection of exclusive decays are 
much easier than those of inclusive ones, although the theoretical understanding of exclusive 
decays is complicated considerably by nonperturbative hadronic form factors. Exclusive 
B — > M£ + £~ decays have been previously studied in ^ [J in the framework of the heavy 
quark effective theory |7|]. Later these decays were also examined for new physics effect 
pi. As is well known, the investigation of the rare decays B — > £ + £~ in future B factories, 
KEK-B, SLAC-B, B-TeV and LHC-B, would provide us of one of the best way to determine 
\V ts \ (and \V td /V td \). 

This paper is organized as follows. In Section 2 we present analytic expressions for 
dT/ds(B — > M£ + £~), M = n, p, K, K*. In Section 3 we study numerical distribution of 
dT/ds(B — > K*£ + £~) and the ratios of branching fractions 

B[B -> p£ + £-)/B(B -> K*£ + £~) and B[B -> tt£ + £-)/B(B -> K£ + £~). 

2 Theory of B -> M£ + £~ (M = n, K, p, AT*) decays 

In the Standard Model the process B — > M£ + £~ (M = tt, if, p, if*) is described at quark 
level by b — > q£ + £~ (q = s, d) transitions and receives contributions from Z and 7 mediated 
penguins and box diagrams.. The QCD corrected Hamiltonian for b — > q£ + £~ decay can be 
written [g 0, [T|] as 



where 14,- are the CKM matrix elements. The operators are given as 



o 1 = 


(qLal^b La ) 




o 2 = 


(qLalftbLp) 




o 3 = 


(qLalfibLa) 


E fafff&fi) 

q'=u,d,s,c,b 


o A = 




E (^VsL) 


o 5 = 


(qLal^b La ) 


E (^a^) 

g'=u,d,s,c,& 
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06 = {qial^Lp) J2 [Qjtp'f^Ra 
q'=u,d,s,c,b 

g 

^7 = 77-£<Za<V (m b R + m q L) baF^ , 
°* = T^q a T^o^{m b R + m q L)b p G a ^ 

107T Z 

9 = ^(q a rLb a )(£^£) , 
O w = ^(q a rLb a )(h^i) , 



where L(R) = | (1 ± 75) are the chiral projection operators. 

Using the effective Hamiltonian in Eq.(l), the resulting QCD corrected matrix element 
for the decays b — > q£ + £~ (q = d, s) can be written as 

M = l=^V*V tb [ci{ f q ltl {\- 75) b£ llx £ + C 10 97m (1-75)^7^ 



'1/ 

- 2C 7 e// g [m, (1 - 75) + m b (1 + Tb )] }> , (2) 



where 

ci{ s = C 9 + (S) . (3) 
In Eq. (2), q v is the four momentum transfer to dileptons, s = q 2 fm\ and 

<% = ^jl + ^ W (S)J + *&(S) . (4) 

The function KJ^ (s) is the one-loop matrix element of Og and Y£ D (s) describes the long 
distance contributions due to the vector J/ip, ift', ■ ■ ■ resonances. The function u (s) repre- 
sents the one gluon correction to the matrix element of the operators Og. Its explicit form 
can be found in |12]|. The explicit forms of the two functions Y^ D (s) and Y^ D (s) are given 
1 as 

Y§ D {s) = g{m c ,8)[ZC x + C % + ZCi + Ct + ZC 5 + Ce\ 

- ^g(l,s)[iC 3 + AC, + 3C 5 + C 6 ] 

- ^(0,s)[C3 + 3C 4 ] + ^[3C 3 + C 4 + 3C 5 + C 6 ] 
" [3Ci + C 2 ] [g (0, s) - (m c , $)] , 



- H-W» cm -W» [X3+Ci+3Ch+c A 

x £ / m ^ +nM V (5) 

The function g (rh q , s) arises from the one loop contributions of the four quark operators 
O x -0 % , i.e., 

8, „ 8 4 2 



g{m q ,s) = --\nm q + — + -y q --(2 + y q )^\l-y q \ 



0(1 - y q ) Mn l±^^ - ^ + 9(y 9 - 1)2 arctan- 



where y 9 = 4m^/s. In Eq. (5), C<® = 3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 . Using m t = 175 GeV, 
nib = 4.8 GeV, m c = 1.4 GeV, a s (m w ) = 0.12 and a s (nib) = 0.22, the numerical values of 
the Wilson coefficients, which we will use in our further numerical analysis, are 

Ct = -0.26, C 2 = 1.11, C 3 = 0.01, 
C A = -0.03, C 5 = -0.03, C 6 = -0.03, 
C 7 = -0.32, C 9 = 4.26, and C 10 = -4.62. 

The values of Cg and C*io are very sensitive to m t . We will neglect the second term in 
Yl D (s), since 3C 3 + C 4 + 3C 5 + C 6 < Cg. Masses, widths and leptonic branching ratios of the 
J p = 1' cc resonances are presented in |T3. Note that as far as short distance effects are 



considered, the w-loop matrix element contribution to the b — > s£ + £~ process is negligible 
due to the smallness of V* s V u b compared to V* b V cs ~ —V t * s V t b, while in the b — > d£ + £~ case, 
the term proportional to V* d V u b is kept. The factor k is chosen to have the value k = 2.3 



1J] to reproduce the rate of decay B — > X s J/ip — > X s £ + £~. The phase of k is fixed, since 
recent experimental data have determined the sign of the ratio of factorization approach 
parameters a 2 /ai and the phase of a± is expected to be near its perturbative value fl5fl . 

At this point, there arises the problem of computing the matrix elements of Eq. (2) 
between the mesons B and M (M = tt, K, p, K*) states. The matrix element (M\Ai\B) 
has been investigated in the framework of different approaches, such as chiral perturbation 



theory fIB| , three point QCD sum rules |I7|| , relativistic quark model fll^H , effective heavy 
quark theory 0, and light cone QCD sum rules [|I~9| j-[[2"2"||. 

The matrix elements for B — > P£ + £~ (P = tt, K) decays can be written in terms of the 
form-factors 

(P{p 2 ) \q^{l - T5 )6| B{ Pl )) = f + (q 2 )( Pl +p 2 ), + f-(q 2 )q» , 



(P(p 2 ) \qior^q v b\ B{p x )) = (p x + p 2 ) M g 2 - (m| - m|)g M 



/r(g 2 ) 

mB + mp 



where q = pi — p 2 - 

The matrix elements for P — > (V — p, K*) decays are defined as follows 



(Vfo.e) |g 7 M(l-7 5 )6|S(Pi)> 



2my 



-— — - + m v) A(g ) + z(pi + p 2 )^ g) 

ffifl + Vfly TUB + Vfly 



+ iqn(e*q)—f- A 3 (q 2 ) - A (q 2 ) 



q 



(7) 



(V(p 2 ,e) \qia^{l + lb )b\B{p l )) = 

4 W e> 2 V 3 T 1 (g 2 ) + 2* \e*{m 2 B - m v ) - ( Pl + p 2 ) ,(e* q)} T 2 (q 2 ) + 



+ 2i{e*q) 



q» - (P1+P2) 



A* 2 

m n — m 



B ~ ""V 



Uq 2 ) 



(9) 



where e is the 4-polarization vector of the \^-meson. Using the equation of motion, the 
form-factor A^{q 2 ) can be written as a linear combination of the form-factors A\(q 2 ) and 
A 2 (q 2 ) (see 0). 

2 mjj + my 2 m B -m v 2 
) = — = Ai(q ) A 2 (q ) , 



2mv 

with the condition A 3 (q 2 = 0) = A (q 2 = 0). 



2my 



Using Eqs. (2), (7), (8), and (9) and summing over the final lepton polarization for 
B -> P£+^ and B -> V7 + £~ decay widths, we get 



c/s 



(p -> Pr 



"B 



2 8 3tt 5 



IW| 2 A 3/2 



x 



2m b C 7 (— C° 9 "f + tf) 



+ 



CW + (<? 2 ) , (10) 



^(P - ^ 
as 



x < 16Am| 



G 2 a 2 m B 
2 12 3tt 5 



l*W| 2 «A 1/a 



|A| 2 +|C<| 2 



+ 2 



-4 [Re (PiP 2 ) + Re (^1^2)] 



|Pi| 2 + |Pi| 
m|A 1 



A + 12rs 



+ 2 



rs 



|P 2 | 2 + |P 2 



miA 2 



rs 



rs 



(11) 



where A = l + r 2 + s 2 — 2r — 2s — 2rs, r = m 2 M /m B , s = q 2 /m 2 B . In Eq. (11) A, Bi, B 2 , C, Di 
and D 2 are defined as 



.4 



C Q e// 



U(g 2 ) 



! ) - — +4^1^) 



mg + my 



q- 



Bl = C e 9 ff (m B + m v )A 1 (q 2 )+AC 7 ^(m 2 B -m 2 v )T 2 (q 2 ) 



B 2 = Cr^^ + AC^(u q 2 ) + ^-^T^ 
tub + rny q z \ m B —m v 



c = (V 



rriB + my 
Di = C w (m B + m v )A 1 (q 2 

A 2 (q 2 ) 



D 2 = C 



10" 



m B + m v 

3 Numerical analysis and discussions 

Now we consider the differential decay widths, dT/dq 2 (B — > (71", p, K, K*) + £ + + The 
hadronic formfactors in the framework of light-cone QCD sum rules have been calculated in 
[r^1-|[22|. For those values of the formfactors, we have used the results of |19| , pOf| , where the 



radiative corrections to the leading twist contribution and SU(3) breaking effects are also 
taken into account. The q 2 dependence of the formfactors can be represented in terms of 
three parameters as 

^ • 

1-ofX- + M4-J 

m B \ m B J 

where the values of parameters F(0), and bp for the relevant decays, B —>■ ir, B —>■ K, 
B — > p and B — > K*, are listed in Table 1. We note that light-cone QCD sum rules method 
is applicable for those decays in the region of ml — q 2 = few GeV 2 , and we found that sum 
rules works very well up to q 2 = 20 GeV 2 . In order to extend our investigation to the full 
physical phase space, we use the above mentioned parametrization in such a way that up to 
q 2 = 20 GeV 2 it successfully reproduces the light-cone QCD sum rules predictions. 

A few words about error analysis in the differential decay rates and their ratios are in or- 
der. In both cases the errors which come from different form-factors are added quadratically 
since they are theoretically independent of each other. Note also that all errors, which come 
from the uncertainties of the b quark mass, the Borel parameter variation, wave functions, 
non-inclusion of higher twists and radiative corrections, are added in quadrature. The un- 
certainty in the ratio is estimated as the half distance between the maximum and minimum 
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F(0) 


cip 


up 


F(0) 




up 






U.ZjU _i_ u.ui 




— n 4i ^ 


n 04 _|_ n nc 


u.uu 




AB- 


^K* 


A?' 


■*p 




Q3 


—0 092 


28 ± 04 


1 18 


281 


A%- 


->K* 


yB- 




34 ± 05 


1.37 


0.315 


46 ± 07 

V 7 • J- V7 — L — V 7 • V J | 


1.55 


0.575 


yB- 


-+K* 


T B- 




0.15 ±0.02 


1.41 


0.361 


0.19 ±0.03 


1.59 


0.615 


T B- 
l 


■*K* 


Tt 


^p 


0.15 ±0.02 


0.28 


-0.500 


0.19 ±0.03 


0.49 


-0.241 


rpB- 

1 2 


■+K* 


Tt 


^p 


0.10 ±0.02 


1.06 


-0.076 


0.13 ±0.02 


1.20 


0.098 


rpB- 


-*K* 


ft 




0.30 ±0.04 


1.35 


0.270 


0.35 ±0.05 


1.37 


0.350 


f?- 


>K 


fr 


>7T 


-0.30 ±0.04 


1.34 


0.260 


-0.39 ±0.05 


1.37 


0.370 




*K 



Table 1: B meson decay form factors in a three-parameter fit, where the radiative corrections 
to the leading twist contribution and SU(3) breaking effects are taken into account (This 



Table is taken from fll9", 20 1). 



values of the ratio. 

In regard to the determination of \Vt s \, we first consider B — > K*£ + £~ . Since the expres- 
sion for B — > V£ + £~ decays contains many form-factors (see Eq. (11)), each with its own 
uncertainty, the error in the differential decay width may be substantial. To reduce these 
uncertainties, it would be better to choose a kinematical region, in which the contributions 
from most of these form-factors will be practically negligible. Therefore, we propose to con- 
sider the end-point region, 16.5 GeV 2 < q 2 < 19.25 GeV 2 . In this region the contributions 
from the terms that are proportional to ~ C-j are much smaller than those from terms pro- 
portional to ~ Cg or ~ Cio- This is due to the fact that the terms which are proportional 
to ~ CV contain 1/q 2 factors. In this region the decay width takes the following form 

G 2 a 2 m\ 



dT 

ds 



2 n 37T 5 



A V2 



X 



Co 



+ \C 



10 



(l±v^)V: 



2 A ± 12rs 



l + <J)(l + A)(l + d) 



(12) 



where 5 represents the contributions from the form-factors V and A 2 . A takes into account 
the magnetic momentum operator contributions (~ Oj), and d(q 2 ) parametrizes the long 
distance effects. In the above-mentioned region, the long distance contributions are taken 
into account in the following way: 



C e J f 



+ \c 



10 1 



Co 



+ \c 



10 1 



1 + d{q 2 
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In order to estimate the uncertainties which arise from d(q 2 ), A(g 2 ) and S(q 2 ), we present 
their q 2 dependences in Figs. 1, 2 and 3, respectively. In Fig. 1 we take into account all 
six J PC = 1 cc resonance states. From this figure we observe that the long distance 
contribution is about ±0.18 of the short distance contribution. As is obvious from Fig. 2, 
the contribution of the magnetic dipole operator is about A ~ — 0.13±0.01, whose behavior is 
observed to be almost independent of q 2 . The same contribution calculated in the framework 
of the HQET gives numbers quite close to our result, in the range —0.18 < A < —0.14 [ O , 
in the same region. In Fig. 3, the dependence of S(q 2 ) on q 2 is depicted. As is clear from 
this figure, the contribution of S to the differential decay width is substantial at lower values 
of q 2 . However, the error due to the uncertainties in the form-factors only causes about 
~ 1% deviation from the case where central values of all form-factors are used. In Fig. 4 
we present the dependence of dT(B — > K*£ + £~)/dq 2 on q 2 and indeed, with the observed 
oscillatory behavior, the uncertainty that d(q 2 ) brings to the B —>■ K*£ + £~ differential decay 
rate is about ~ 8%, although its leading uncertainty to the perturbative contribution is 
about ±0.18. 

We note that the extraction of |Vtd/Ks| from the decay widths B — > p£ + £~ and B — > 
K*£ + £~ in the low invariant mass region, for example 2 GeV 2 < q 2 < rn 2 ^, becomes more 
problematic. In this region the contribution of the magnetic dipole operator O7 is large. 
Therefore the theoretical predictions of the decay widths of the processes B — > p£ + £~ and 
B — > K*£ + £~ have significant uncertainties in this region, since the form-factors 7\, T2 
and T3 play an essential role. In our opinion, it is better to investigate the B — > tt£ + £~ and 
B — > K£ + £~ decays in this low momentum region for the determination of the ratio |K<i/Vt s |, 
since both decays are described only by the two form-factors /+ and fx (when lepton masses 
are neglected). Various theoretical schemes predict that /+(0) ~ — /t(0) (see ||l3 _ lHl)- m 



this region of q , that is far from q^axi ^ e B* pole contribution |L4| is unlikely to violate 
the SU(3) symmetry predictions. 

We now show the possibility for extracting |T4d/T4s| from the ratio 

d 

dq 2 



dR 

dq 2 



B(B -> p£ + £-)/B(B -> K*rr) 



(13) 



and the ratio 
dR x 



d 



dq 2 dq 2 



B(B 



~)/B(B -> K£ + t 



td 



(14) 



Performing the integrations over q 2 , from 16.5 GeV 2 to 19.25 GeV 2 in Eq. (13) and from 2 
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GeV 2 to 9 GeV 2 in Eq. (14), i.e. 



R(16.5 < q 2 < 19.25 GeV 2 ) - B ( B ^ P £ 



Ri(2 < q 2 < 9 GeV 2 ) 
we get the numerical results 



g(g -» ttI+1 



R = 0.85 ±0.06 , 
Ri = 0.77 ±0.06, 

(normalized to IV^/V^I 2 ). 

At this point we estimate the number of expected events of the types B — > 
and -B — > P£ + £~ in experiments at future S factories. Future symmetric and asymmetric 
B factories of electronic and hadronic colliders should produce much more than 10 9 B B 
mesons per year by the year 2010. Assuming 10 9 B mesons effectively reconstructed, the 
number of expected events in the corresponding kinematical regions (16.5 GeV 2 < q 2 < 
19.25 GeV 2 for B -> V£ + £~ decay and 2 GeV 2 < q 2 < 9 GeV 2 for B -> P£ + £~ decay) are 
calculated to be 



N(B -» 


K*£ + £~ 


) = 


-> 


r) x io 9 ~ 140 


N(S - 


-> #£ + r 


) = 


^«+^ 


") x 10 9 ~ 735, 


N(S 




) = B(S- 




■) x 10 9 ~ 5, 


N(S - 




) = B(S- 




") x 10 9 ~ 28. 



As can be seen easily, for determination of \V t d\ or |V^/Vt s | we need much more than 10 9 B 
mesons produced, and it will be only possible at future hadronic B factories, such as B-TeV 
and LHC-B, where the double lepton triggering helps high reconstruction efficiencies with 
more than 10 11 B mesons per year produced. From the above results we conclude that these 
decays have a good chance to be detected at future B factories. 

The problem of accurate determination of CKM matrix elements with different methods 
receives special attention in connection with the fact that the next generation of B meson 
decay experiments will be a test of the flavor sector of the SM at high precision as well 
as allowing the determination of V t d and with very high accuracy. Simultaneous deter- 



minations of CKM angles and phases [ffl, 23 1 would be extremely important to check the 



consistence within the SM and to search for the hints of new physics beyond the SM. It is 
well known that the experimental investigation and detection of exclusive decays is much 
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easier than those of inclusive ones, while the theoretical understanding of exclusive decays is 
much more complicated due to nonperturbative hadronic form factors. We have found that 
for the determinations of \V ts \ (and |Vtd/Vt s |), due to the character of the relevant hadronic 
form factors, the kinematic region of small q 2 is more useful for the B — > P£ + £~ decays, 
while the kinematical region of q 2 near the end-point is more suitable for the B — > Vi + i~ 
decays. And we estimated numerically \V ts \ and the ratio |V^/VJ S | from the ratios of those 
decays with a theoretical uncertainty of ~ 10% . 
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Figure Captions 

1. Dependence of d on q 2 for B — > K*£ + £~ decay. Dash-dotted d(q 2 ) = line corresponds 
to the short distance contributions. 

2. Dependence of A on q 2 for B — > K*£ + £~ decay. Dotted and dash-dotted curves cor- 
respond to the cases when the uncertainty is added and subtracted from the central values 
of all form-factors, respectively. 

3. Dependence of 5 on q 2 for B — > K*£ + £~ decay. Dotted and dash-dotted curves cor- 
respond to the cases when the uncertainty is added and subtracted from the central values 
of all form-factors, respectively. 

4. Dependence of dT{B — > K*£ + £~)/dq 2 on q 2 in units of 

(l + v^) 2 • 

In this figure the wavy curve numbered as 2 represents the case that takes into account all 
1 states of cc resonances for the central values of all form-factors, while the wavy curves 
numbered as 1 and 3, correspond to the cases when the uncertainty is added and subtracted 
from the central values of all form-factors, respectively. The dash-dotted line represents the 
contribution of only the three lightest resonances. The dotted curve, on the other hand, is 
for the perturbative result, i.e., d(q 2 ) = 0. 



^ g m B | V T/ *l 2 
2 ll 37r 5 



Cc 
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10| 



12 



References 

[1] For example, for the relation between V td and V ub , see T. M. Aliev and C. S. Kim, Phys. 
Rev. D58 (1998) 013003. 

[2] M. Gronau and J. L. Rosner, Phys. Lett. B376 (1996) 205. 

[3] G. Buchalla and A. Buras, Phys. Rev. D54 (1996) 6782. 

[4] C. S. Kim, T. Morozumi and A. I. Sanda, Phys. Rev. D56 (1997) 7240. 

[5] W. Roberts, Phys. Rev. D54 (1996) 863. 

[6] W. Roberts and F. Ledroit, Phys. Rev. D53 (1996) 3643. 

[7] C. Greub, A. Ioannissian and D. Wyler, Phys. Lett. B346 (1995) 149. 

[8] N. Isgur, M. B. Wise, Phys. Rev. D41 (1990) 151; 
M. B. Wise, Phys. Rev. D45 (1992) 2188; 

G. Burdman and J. F. Donoghue, Phys. Lett. B280 (1992) 287; 
L. Wolfenstein, Phys. Lett. B291 (1992) 177; 
T. M. Yan et. al. Phys. Rev. D46 (1992) 1148; 

G. Burdman, Z. Ligeti, M. Neubert and Y. Nir, Phys. Rev. D49 (1994) 2331. 

[9] B. Grinstein, M. J. Savage and M. B. Wise, Nucl. Phys. B319 (1989) 271. 

[10] M. Misiak, Nucl. Phys. B398 (1993) 23; Erratum, ibid. B439 (1995) 461. 

[11] M. Ciuchini, E. Franco, G. Martinelli, L. Reina and L. Silvestrini, 
Phys. Lett. B 316 (1993) 127; 

M. Ciuchini, E. Franco, G. Martinelli and L. Reina, Nucl. Phys. B 415 (1994) 403; 
G. Cella, G. Curci, R. Ricciardi and A. Vicere, Nucl. Phys. B 421 (1994) 41; 
ibid. Phys. Lett. B 325 (1994) 227. 

[12] A. J. Buras and M. Miinz, Phys. Rev. D 52 (1995) 186. 

[13] Particle Data Group, Phys. Rev. D54 (1996) I. 

[14] Z. Ligeti and M. Wise, Phys. Rev. D53 (1996) 4937. 

[15] T. E. Browder, K. Honscheid and D. Pedrini, 
Ann. Rev. Nucl. Part. Sci. 46 (1996) 395. 



13 



[16] R. Casalbuoni et al, Phys. Rep. 281 (1997) 145. 

[17] P. Colangelo, F. De Fazio, P. Santorelli and E. Scrimieri, 
Phys. Rev. D53 (1996) 3672; Phys. Lett. B395 (1997) 339. 

[18] W. Jaus and D. Wyler, Phys. Rev. D41 (1990) 3405; 

D. Melikhov, N. Nikitin and S. Simula, |hep-ph/ 9 704268] (April 1997). 



[19] P. Ball, |hep-ph/9802394| (February 1998). 



[20] P. Ball, V. M. Braun, |hep-ph/ 980542^ (May 1998). 



[21] T. M. Aliev, A. Ozpineci and M. Savci and H. Koru, Phys. Lett. B400 (1997) 194. 

[22] T. M. Aliev, A. Ozpineci and M. Savci, Phys. Rev. D56 (1997) 4260. 

[23] A. S. Dighe, M. Gronau and J. L. Rosner, Phys. Rev. D54 (1996) 3309; C. S. Kim, D. 
London and T. Yoshikawa, Phys. Rev. D57 (1998) 4010. 



14 



0.2 



0.1 



0.0 - 



-0.1 



-0.2 



-0.3 




q 2 (GeV 2 ) 



Figure 1: 



-0.150 




-0.170 - 



15.0 



16.0 17.0 

q 2 {GeV 2 ) 



18.0 



19.0 



Figure 2: 



15 



0.05 - 



I 0.00 

+ 



-0.05 



-0.10 




15.0 



16.0 17.0 

q 2 {GeV 2 ) 



18.0 



19.0 



Figure 3: 



0.60 



+ 

t 

r I ^ 
^3 | "T3 



0.45 - 



0.30 



0.15 



0.00 




15.0 



16.0 17.0 

q 2 {GeV 2 ) 



18.0 



19.0 



Figure 4: 



16 



